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The magnetic structure of multiferroic M3V2O8 has been investigated using non-resonant X-ray 
magnetic scattering. Incident circularly polarized X-rays combined with full polarization analysis of 
the scattered beam is shown to yield high sensitivity to the components of the cycloidal magnetic 
order, including their relative phases. New information on the magnetic structure in the ferroelectric 
phase is obtained, where it is found that the magnetic moments on the "cross-tie" sites are quenched 
relative to those on the "spine" sites. This implies that the onset of ferroelectricity is associated 
mainly with spine site magnetic order. We also demonstrate that our technique enables the imaging 
of multiferroic domains through polarization enhanced topography. This approach is used to image 
the domains as the sample is cycled by an electric field through its hysteresis loop, revealing the 
gradual switching of domains without nucleation. 

PACS numbers: 75.25.-j, 75.85.+t, 77.80.Dj 



I. INTRODUCTION 

X-ray imaging continues to find growing utility driven 
in part by the availability of intense beams from syn- 
chrotron and free-electron sources, and in part by a grow- 
ing awareness of how to exploit various contrast mech- 
anisms and imaging modes. Imaging techniques range 
from those that depend on monitoring the absorption, 
or phase shift 1 3 , to diffraction of an x-ray beam, in- 
cluding the possibilities offered by coherent beams 4 ^ and 
holograph}*^. Here we introduce a technique where 
imaging proceeds by measuring the polarization state of 
a magnetically scattered photon through determination 
of its Stokes parameters. We demonstrate that when 
applied to multiferroic M3V2O8 it provides topographic 
images of the spatial distribution of magnetic domains 
as the ferroelectric domains are cycled through a hys- 
teresis loop by an applied electric field. Our technique 
has potential for imaging domains in multiferroic devices 
and other classes of correlated electron systems which are 
characterised by unconventional or coupled order param- 
eters. 

The discovery of strongly coupled magnetic and ferro- 
electric order parameters in TbMnO^, and subsequently 
in other families of compounds^^, has led to a renais- 
sance of interest in magnetoelectric materials. It is now 
established that TbMn03 typifies a new class of mate- 
rials in which the onset of ferroelectricity is driven by 
the formation of non-collinear, cycloidal magnetic order 
which breaks spatial inversion symmetry. One of the out- 
standing challenges in this field is the development of 
methods capable of imaging multiferroic domains, and 



most especially their evolution under applied external 
fields. However, the imaging of ferroelectric and mag- 
netic domains has largely remained two separate fields, 
with a few exception s 11 ! 12 . Imaging of ferroelectric do- 
mains is well established, and can be achieved through a 
variety of probes including X-ray charge scatterin g 13 ! 14 , 
and atomic force microscop}^. In contrast, the imag- 
ing of antiferromagnetic domains has emerged in more 
recent times, driven by the availability of highly brilliant 
X-ray beams from synchrotron sources. The high spa- 
tial resolution attainable with such sources has enabled 
new imaging methods based on either the absorption 1 ^ or 
scattering ^ 17 i 18 of an X-ray beam through various pro- 
cesses which yield sensitivity to the antiferromagnetic or- 
der. 

For the new class of multiferroics exemplified by 
TbMnOs, the potential of X-rays for providing infor- 
mation on the multiferroic state has been highlighted 
by several studies^—. In particular, we have demon- 
strated earlier that the handedness of circularly polarised 
X-rays naturally couples to the handedness of cycloidal 
spin structures, and that the change in polarization of 
an X-ray photon undergoing a non-resonant magnetic 
scattering process encodes detailed information on the 
cycloidal magnetic order—. The polarization state of 
a photon may be described by the Stokes parameters 



P=(Pi,P 2 ,P 3 



\27,28 



, where a linear polarization analyser 



may be used to determine the Stokes parameters P\ and 
P2. By definition P\ describes the component of the lin- 
ear polarisation in and perpendicular to the scattering 
plane, P2 the oblique linear polarisation and P3 the cir- 
cular polarisation. Here we demonstrate the imaging of 
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FIG. 1: Crystallographic and low temperature incommensurate (LTI) magnetic structure of the nickel sublattice in Ni3V20s. 
(a) Schematic of the two nickel sites: spine (red) and cross-tie (blue) (b) Projection onto the a-b plane of the proposed LTI 
magnetic structure of the spine sites m SP = (1.6 =b 0.1, 1.3 =b 0.1, 0)/xb and the cross-tie sites m CT = (2.2 ± 0.1, 1.4 ± 0.1, 0)/xb, 
as deduced from neutron diffraction measurements^. 



magnetic domains, and hence the ferroelectric domains, 
in multiferroic M3V2O8 where the imaging contrast de- 
rives from the use of circularly polarized X-rays to ex- 
cite non-resonant magnetic scattering processes, and to- 
pographic imaging is performed by mapping the spatial 
variation of the Stokes parameters. 



M3V2O8 (NVO) is an example of the new class of mag- 
netoelectric multiferroics. It is a magnetic insulator with 
a structure characterised by planes of Ni 2+ S = 1 spins 
arranged on a buckled kagome staircase structure, result- 
ing in two inequivalent Ni sites: "cross- tie" and "spine" - 
type, as shown in Fig. [2(a). Whilst the ideal kagome lat- 
tice is the canonical example of a frustrated system, the 
deviation from this geometry in NVO leads to additional 
interactions relieving the frustration and producing a rich 
magnetic phase diagra m 19 i 29 i 3Q . Below Tl = 6.3 K, in the 
low temperature incommensurate (LTI) phase, the spins 
on both site types are reported to be arranged in a cycloid 
in the a — b plane with propagation vector (r w 0.27 0), 
breaking spatial inversion symmetry. It is in this phase 
that a spontaneous electric polarisation develops parallel 
to the 6-axis^. The magnetic structure in the LTI phase 
derived from the neutron diffraction studie a 19 i 31 is shown 
in Fig. [2(b). It should be noted that in these studies it 
was not possible to uniquely determine the phase rela- 
tionship between the two sites. Recent polarised neutron 
diffraction measurements performed in an applied electric 
field E have demonstrated that the handedness of the cy- 
cloidal order can be switched by reversing although 
no information on the spatial distribution of domains was 
obtained. 



In this work we present the determination of the mag- 
netic structure of M3V2O8, refining the spin moments on 
the two nickel sites and their phase relationship, and a 
demonstration of the imaging of the magnetic domains 
and their electric field control. 
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FIG. 2: Schematic of the experimental set up, where for a 
given magnetic satellite the Stokes parameters are determined 
by measuring the scattered intensity as a function of the angle 
77. 



II. EXPERIMENTAL DETAILS 

NVO is orthorhombic with high temperature crystal- 
lographic structure Cmce (number 64 in the Interna- 
tional Crystallographic Tables) and lattice parameters 
a = 5.92 A, b = 11.38 A, and c = 8.24 A. The single 
crystal used for this experiment was grown from a BaO- 
V2O5 flux^S, and cut and mounted to give a specular 
a- face, with surface dimensions 600 /imx900 fim. Exper- 
iments were performed on the ID20 Magnetic Scattering 
beamline 32 at the European Synchrotron Radiation Fa- 
cility in Grenoble, using a monochromatised X-ray beam 
at an energy of 7.45 keV, considerably below the nickel 
if -edge to avoid interference with the resonant signal and 
to enable the separate determination of the spin and or- 
bital components of the magnetic moments. The experi- 
mental setup is presented in Figure [2j 

A polarisation analyser, assembled on the detector arm 
of a six circle diffractometer, was used to discriminate 
the scattered photon polarisation by using a high quality 



3 



Au(222) analyzer with mosaic spread of ~ 0.22°. At this 
energy the cross-talk between the a' (rotated, 77 = 0°) 
and tt' (unrotated, r] = 90°) channels was approximately 
0.3%. This enables the extraction of the Stokes param- 
eters by fitting the dependence of the scattered light 
/ = Iq(1 + Pi cos(2^) + P2sin(2^)) on the angle 77 of 
the analyser assembly about the scattering vector k 27 i 28 . 

An in- vacuum phase plate setup was used in quarter- 
wave mode to convert the incoming linearly polarised 
light (tt) into the circular left (LCP) and right (RCP) 
components, by using a high quality (110) diamond sin- 
gle crystal plate (1.2 mm thick) in Laue geometry, and 
oriented at about 45° with respect the incident polarisa- 
tion 7T, in order to have the (111) reflecting plane normal 
to the surface 33 . The Stokes parameters of the incident 
beam were carefully checked before and after every set of 
scans by using the same polarisation analyser, in order 
to have a good control of the incident circular photon 
polarisation (in this case |Pi| « IP2I ~ 0.008(4), hence 
circular beam polarisation close to 99%, assuming no de- 
polarisation of the beam such that P3 = \f(X — P\— P| ))? 
which is very sensitive to the beam position stability. 

Measurements were performed by recording rocking 
curves of the analyser at a series of values of 77. The in- 
tensity I(r]) was then obtained by numerically integrating 
the individual peak shapes for LCP and RCP incident. 
To correct our data for the diffuse Thomson scattering 
background, arising from the sigma component of the cir- 
cular polarisation, measurements were also performed off 
the magnetic reflection (AO = 0.2°) and then the back- 
ground I(r]) was subtracted from the signal at the reflec- 
tion. 

An electric field of up to 2 kV/mm was applied by 
means of two copper electrodes glued by high conduc- 
tive silver paste on the sample 6-surfaces. The sample 
stick was inserted into an orange cryostat allowing us to 
access a temperature range of 2 — 300 K. NVO is an insu- 
lator, and in spite of the use of 4 He exchange gas, given 
an incident photon flux of ~ 10 13 ph/sec, it was neces- 
sary to attenuate the incident X-ray beam to reduce the 
transmission by a factor of five, in order to maintain the 
sample within the LTI phase. 
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FIG. 3: Temperature dependence of the scattering intensity 
(upper) and wave-vector position (lower panel) of the mag- 
netic satellite (5 — r,l,0) measured in the ty — a' channel. The 
lines^ define the boundaries of the low temperature incom- 
mensurate phase of interest. The inset in the top panel shows 
a typical reciprocal space scan along the /i-direction of the 
same satellite, measured at T = 4 K. 



interest was clearly defined, and our measurements were 
performed at an optimised temperature, taking great 
care to avoid thermal instabilities which could lead to 
crossing the phase boundaries. 



A. Magnetic Structure Refinement 



III. EXPERIMENTAL RESULTS AND 
ANALYSIS 

Figure [3] shows the temperature dependence of the 
magnetic reflection (5 — r,l,0) obtained from scans mea- 
sured in the rotated tt — a channel along the reciprocal 
space /^-direction. The sharp change seen in the inten- 
sity of the scattering at T = 3.2 K, obtained by numeri- 
cal integration, indicates the transition between the low 
temperature commensurate (LTC) and incommensurate 
(LTI) phases, whilst the discontinuity seen in the evolu- 
tion of the propagation wave- vector at T = 6.3 K occurs 
at the transition between the high (HTI) and low temper- 
ature incommensurate phases. Hence the LTI phase of 



Our investigation of the magnetic structure was facil- 
itated by the fact that an applied electric field can cre- 
ate a magnetic domain with a single handedness through 
the magneto-electric coupling^, and was accomplished 
by performing a full Stokes analysis of the scattered in- 
tensity, a technique which, as established for TbMnOa^, 
provides information on the magnetic structure that is 
not accessible to neutron diffraction experiments. The 
field cooling procedure for performing the measurements 
was to apply —1.4 kV/mm across the sample as it was 
cooled from 20 K, in the paramagnetic phase, to 4.2 K, 
in the LTI phase. The electric field was then removed 
before illuminating the sample with X-rays. 

Figure |4] shows the Stokes dependence measured for 
three different reflections: (5 + r, 1, 1), (5 + r, 1, 0) 
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FIG. 4: The Stokes dependence of the scattered intensity 
from magnetic satellites on analyser rotation angle r\ was de- 
termined by the data measured in the LTI phase with left 
(LCP) and right circularly polarised (RCP) X-rays after field 
cooling with E = 1 kV/mm along —b. The solid lines repre- 
sent calculations with spin moments of raf p =1.6, raf p =1.3, 
rria T =0, and raf T =0, in units of /is, and a domain popu- 
lation of 87%. The dashed lines correspond to the neutron 
mode^ with finite cross-tie moments (Fig.[TJb)), with an as- 
sumed phase difference of ty/2 between the spine and cross-tie 
moments, and where the intensities are scaled to that of the 
upper panel. 



and (5 + r, -1, 0) with both LCP and RCP X-rays in- 
cident, as red and green markers respectively. Attempts 
to observe a signal at a fourth reflection: (4 + r, 1, 0) 
were unsuccessful. In all three panels one notices that 
the maxima and minima in the integrated intensities lie 
close to 45° and 135°, where their positions are reversed 
for LCP vs RCP incident. In addition the large varia- 
tion in the intensities as a function of rj indicates that 
the X-rays have been effectively converted on scattering 
from being circularly to linearly (in this case obliquely) 
polarised, i.e. P2 > P\ > P3. This differs from the 
behaviour seen for TbMnOa, in which the scattered X- 
rays retained a partially circular polarisation stated. For 
(5 + t, 1, 1) in the first panel, the maximum intensity for 
LCP is greater than that for RCP, whilst in the lower two 



panels the maximum intensities are approximately equal 
for the two circular polarisations. 

A particular feature of NVO is that by taking advan- 
tage of the structure factor one may selectively probe the 
spine and cross-tie moments independently or in combi- 
nation. Magnetic reflections can be classified according 
to their indices (hkl): type 1 - h odd, k odd, I odd - 
cross-tie moments cancel such that only the spine site 
moments contribute, type 2 - h even, k even, I odd - 
the reverse is true, such that the cross-tie moments are 
singled out, and type 3 - h odd, k odd, I even - both 
moment types contribute to the scattering. The three 
different peak type scattering amplitudes are: 



fi oc 8cos(27r&A^)(mf p - vyam% 
h oc 4(m^T - i 7 amf T ) • B, 
h oc (8/3sin(2^A| p )(mf p -i 7 amf p ) 
+4exp(i</>)(m^ T - i7«mf T )) • B, 



5P )-B, 



(1) 



where a = ±1 gives the sign of r in the propagation 
wavevector, f3 = (— l) n +^ for h = 2n + 1 and I = 2p, 
7 = ±1 for the two different domains, Aj? p = 0.13 is the 
fractional coordinate along the 6-axis of the Wyckoff po- 
sition occupied by the spine ions, describing the buckling 
of the Kagome planes, and the vector B is that of Blume 
and Gibbs^i, containing the polarisation state of the inci- 
dent and scattered X-rays for the spin contribution. Thus 
the different magnetic reflections we investigated can be 
identified as : type 1 -(5 + r, 1, 1), type 2 - (4 + r, 
1, 0) and type 3 - (5 + r, +(— )1, 0), where the sign of 
+ (— )1 enables us to extract information regarding the 
phase relationship <fi between the two moment types. 

Starting with the pure spine- type only reflection (5 + r, 
1, 1), we compare our data for both LCP and RCP with 
a calculation for the magnetic structure proposed in an 
earlier neutron diffraction study 19 , with total magnetic 
moments: 



m 



SP 



m 



SP 



1.6 /is, 
1.3 fl B - 



Our calculations use these total magnetic moments as 
spin moments only on the Ni ions, having assumed that 
the orbital moment is negligible^, and find an excellent 
agreement. Our fit, which is shown in the upper panel of 
Fig. [4] as a solid line, allowed us to extract the domain 
population which we found to be 87(4)% domain type 1, 
corresponding to the moments rotating clockwise in the 
a — b plane as we move along the a axis, when observed 
from -fc. 

The absence of an observed signal at the cross-tie site 
only reflection suggests that any moment on this site is 
very small. For the third type of reflection, exemplified 
by the twin reflections (5 + r, ±1,0), to aid the under- 
standing of the data and to investigate the cross-tie site 
moment and the phase relationship between the two mo- 
ment sites, we can make use of a simplifying model. 

We stress that the purpose of this simplified model is 
to allow us to understand the qualitative behaviour of 
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the data shown in Fig. [4] only. The starting assumptions 
for the simplified model are: 1) a single domain, 2) an 
analyser Bragg angle of 6 = 45°, 3) the scattering plane is 
perfectly horizontal, and 4) within the scattering plane, 
the angle between the 6-axis and the scattering vector is 
45°, the scattering amplitudes in the a' and tc' channels 
can be written as: 

U> oc [kS a + C b sin - C a cos 0] 

+i[kSb — Cb cos <j) — C a sin <j)} , 

fir' oc [kSb + Cb cos (j) + C a sin </>] 

+ie[kS a — Cb sin — C a cos </>] , 



where 



S a 

s b 

C a 



SP\^SP 



8sin(27rA^)m 
8sin(27rAf p )raf p , 



4m 



CT 



C 6 = 4mf T , 

fc = ±1 differentiates between the two reflections (5 + 
r, ±1, 0), and e = ±1 refers to the handedness of the in- 
cident light. The average intensity for the two reflections 
can be expressed as: 



\f*>\ 2 + \f*> 



\kS a + Cosine 
| k Sb — Cb cose 
\kSb + Cb cose 
\kS a - C b sine 



C a COS ( 

C a sin <; 
C a sin <; 
C a cos <j>\*. (2) 



It is noticeable that e, the sign of the incident polarisa- 
tion, is now absent, consistent with the average intensity 
for LCP and RCP being the same. This arises from the 
spin moments lying in the scattering plane, in contrast 
to when measurements are made of the spine-type only 
reflection. If the observational constraint of equivalence 
between ±fc, i.e. between panels two and three of Fig. |4j 
is applied to Eq. (j2]) this requires that C a cos </> = and 
Cb cos (j) = and hence a phase difference of <p = ir/2. 
This results in further simplifying the scattering ampli- 
tudes to: 



U> cx [kS a + C b ]+i[kS b -C a ], 
U> cx [kS b + C a ]+ie[kS a -C b ]. 



(3) 



If one then considers the observation that the scatter- 
ing is predominantly described by P2, with Pi close to 
zero, where 



Pi 
P2 



\fa 



\h^ 2 + \h^ 2, 
\h' + h'\ 2 -\h 



(4) 



this imposes additional constraints. Firstly for Pi « 
this requires that \f&'\ ~ \f%' |, and hence that S a ~ Sb 
and C a ~ Cb, which is broadly consistent with previously 
published values. Meanwhile P2 is not equal to zero for 
\f#, + fft,\ zfz \ f~, — f~,\^ expressions for which are: 



e = +1 : \U> + h>? = \(S a + S b ) + (C a + C b )\ 2 + \(S a + S b ) - (C a + C b )\\ 

\f*> ~ U>\ 2 = \(Sa ~ S b ) + (C a - C b )\ 2 + \{S a ~ S b ) - (C a - C b )\ 2 , 

e = -1 : \U> + h'? = \(S a ~ S b ) + (C a - C b )\ 2 + \(S a - S b ) - (C a - C b )\ 2 , 

\U> - U'\ 2 = \(Sa + S b ) + (C a + C b )\ 2 + \(S a + S b ) ~ (C a + C b )\ 2 , (5) 



i.e. they are independent of the sign of k. Equation (j5j) 
clearly shows that reversing the handedness of the circu- 
lar polarisation e results in switching the sign of P2 . Since 
we have already established from Pi w that S a ~ Sb 
and C a ~ Cb, then (S a — Sb) and (C a — Cb) are close to 
zero, whilst \(S a + Sb) + (C a + C )\ is strong, such that 
there is a high contrast between \f&> + fa \ and \ — |, 
resulting in large values for P2, which is reversible de- 
pending on e, consistent with the key observations from 
the data shown in Figure [H 

Using the moment magnitudes and intensity scaling 
factor extracted for the spine- type sites given above, a 
fit was made simultaneously to the different data for the 
twin reflections using the full model given in Eq. (pQ), 
including the full six-circle geometry, and this is shown 
by the solid lines in the lower two panels of Figure 2J 



The scattered signal is found to be very sensitive to the 
value of the phase difference between the moments on the 
spine and cross-tie sites, with the similarity in the scat- 
tered intensities for the two reflections implying a phase 
difference of 0.50(2)7r. Calculations for different mag- 
netic structure models allow us to provide estimates of 
m^ T = 0.55(50) and mf T = 0.35(35) (units of fj, B ), with 
a 87(4) % domain type 1 population, where the large er- 
ror bars reflect the relative insensitivity to the cross-tie 
moments. Nonetheless, the data and model calculations 
displayed in Fig. |4] underline the remarkable sensitivity 
of the technique to the magnetic amplitudes and phase 
relationships of complex magnetic structures. It is note- 
worthy that the observed cross-tie moments are consid- 
erably smaller than the published values obtained from 
neutron diffraction^ (calculations for which are shown in 
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FIG. 5: (a) Image of the sample sandwiched between the Cu electrodes used to form a capacitor. The dotted mesh superimposed 
on the image indicates points at which the domain population was determined. The ellipse indicates the size of the X-ray beam, 
(b) Variation of the Stokes parameters as a function of domain volume fraction, demonstrating the sensitivity to P2 (blue line) 
and insensitivity to Pi (red line), allowing a contrast to be measured between the two cycloidal domains, (c) Images of the 
domain populations as a function of applied electric field for the magnetic reflection (5 + r, 1, 0), where the colour represents 
the percentage of the clockwise and ant i- clockwise magnetic cycloidal domains. 



Fig. 2] as dashed lines), although the analysis of more re- 
cent neutron diffraction measurements seems consistent 
with the possibility of there being no moment on the 
cross-tie nickel ions [I. Cabrera et al, unpublished]. Thus 
we may conclude that the cross-tie moments are most 
likely negligible, from which it follows that spatial inver- 
sion symmetry is broken by the non-collinear magnetic 
structure of the spine-sites alone. 



B. Domain Imaging 

The data in Fig. [4] also suggest a novel route to imag- 
ing multiferroic domains. The reversal evident in the sign 
of the Stokes parameter P2 when switching the incident 
X-ray polarization, would also occur for fixed incident 
polarization if the cycloidal magnetic domain switched 
its handedness^. It follows that the domain population 
can therefore be determined at different points in the 
sample by determining the value of P2 (see Fig. [5fb)). 
(Pi and P3 are relatively insensitive to domain volume 
fraction.) For the imaging the incident beam spot was 
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FIG. 6: Position dependent magnetic domain population hys- 
teresis loops obtained from the domain images shown in Fig. 
[5] averaged over different areas of interest: (a) entire sample 
area; (b) edge; (c) central region. In (c) the data are com- 
pared with bulk measurements of the electric polarisation 31 
where an appropriate scaling has been applied. 



reduced to 0.25 x 0.38 mm 2 to maximise the resolution 
whilst still giving a reasonable count rate in the detec- 
tor (Fig. [5ja)). The sample was then rastered through 
the beam and at each point, the scattered intensity was 
recorded for the analyser setting angles of r] = 45° and 
135°. The Stokes parameter P2 was evaluated from these 
measurements, hence providing the percentage of left or 
right handed magnetic cycloidal domains. To establish 
the evolution of the domains as a function of electric field 
a different field application protocol was employed. Ini- 
tially the sample was cooled to 5 K with no voltage across 
the sample. Then between each set of measurements, 



whilst remaining at 5 K, the voltage across the sample 
was ramped up to the required level. Images of the do- 
mains are shown in Fig. [5fc) as a function of electric 
field applied, demonstrating clearly that we were able to 
resolve inhomogeneities in the domain populations with 
good resolution. 

Inspection of the individual images reveals that for zero 
applied electric field the distribution of the two domains 
is not equal, but instead shows a preference towards do- 
main 2, i.e. moments rotating anti-clockwise. This may 
either be a surface strain effect or some memory effect of 
an earlier domain state induced by field cooling in pre- 
vious experiments, as has been observed in MnWC^ 35 . 
Applying an increasing electric field along —b leads to a 
reversal of the domain populations, with the production 
of almost a mono-domain 1 state for E = +270 V/mm. 
The evolution of the domain populations is gradual, with 
the boundary between the two shifting as a function of 
applied electric field, as opposed to the nucleation of do- 
main 1 within domain 2 leading to a more randomised 
distribution. When the electric field direction is then re- 
versed the domain distributions are little changed until 
between —250 V/mm and —370 V/mm the distribution 
is again reversed with a preference towards domain 2. 
Whilst the domain populations are preserved on return- 
ing to E = V/mm, due to time constraints no ad- 
ditional data were collected for applied positive electric 
fields. 

One way to quantitatively render the domain homo- 
geneities is to average the images over different areas of 
interest to produce position dependent magnetic domain 
population hysteresis loops. Figure [6fa) shows the do- 
main population averaged over the entire sample. On 
average the excess population of domain 2 in zero field is 
quickly reversed in a positive applied field, with domain 
1 prevailing such that even large negative fields are in- 
capable of restoring an excess of domain 2. In Fig. [6jb), 
the response from the bottom left-hand corner is plot- 
ted, where it is clear that close to the edge of the sample 
there is strong pinning of domain 1 with it dominating 
at all fields. Imaging allows these edge effects to be ex- 
cluded from the data, (Fig. Efc)), and restores more of 
a symmetrical character to the electric field dependence 
of the magnetic domain population. This allows us to 
compare our data with the electric polarisation measure- 
ments of Cabrera et al^i, as shown in Fig. [6jc), where 
we find excellent agreement. This establishes the link 
between the magnetic domains and the electric polarisa- 
tion in multiferroic NVO, proving that by imaging the 
magnetic domains we are in effect also imaging the fer- 
roelectric domains. 



IV. CONCLUSION 

In conclusion, we have refined the magnetic structure 
of M3V2O8, confirming the cycloidal ordering of the Ni 
spins on the spine sites whilst finding that there is no 
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ordered moment on the cross-tie sites, indicating that 
it is the magnetic order on the spine sites alone that 
breaks inversion symmetry. Further, we have successfully 
demonstrated how polarisation enhanced X-ray imaging 
(PEXI), through the contrast provided by the magnetic 
structure, enables the imaging of the magnetic domains 
in NI3V2O8 which endow this material with its multifer- 
roic properties. 

Comparison of the magnetic domain population hys- 
teresis loops with bulk electric polarisation measure- 
ments reveals the coupling between the magnetic and 
ferroelectric domains. This opens the prospect of using 



this technique to image multiferroic domains in related 
materials, and even in operational devices. With the im- 
proved focussing offered by beamlines currently under 
development the spatial resolution for this technique will 
be reduced to 100 nm and below, at which point it will be 
possible not only to image the domains themselves, but 
also the structure of the domain walls, the engineering 
of which ultimately determines the usefulness of a given 
material. 

Many thanks to everyone who has helped us with the 
experiments, especially A. Fondacaro, J. Herrero-Martin, 
C. Mazzoli, V. Scagnoli and A. C. Walters. 
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